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Abstract. We briefly review some equilibrium and nonequilibrium properties of systems with long-range 
interactions. Such systems, which are characterized by a potential that weakly decays at large distances, 
have striking properties at equilibrium, like negative specific heat in the microcanonical ensemble, tem- 
perature jumps at first order phase transitions, broken ergodicity. Here, we mainly restrict our analysis to 
mean-field models, where particles globally interact with the same strength. We show that relaxation to 
equilibrium proceeds through quasi-stationary states whose duration increases with system size. We pro- 
pose a theoretical explanation, based on Lynden-BelPs entropy, of this intriguing relaxation process. This 
allows to address problems related to nonequilibrium using an extension of standard equilibrium statistical 
mechanics. We discuss in some detail the example of the dynamics of the free electron laser, where the 
existence and features of quasi-stationary states is likely to be tested experimentally in the future. We 
conclude with some perspectives to study open problems and to find applications of these ideas to dipolar 
media. 
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1 Introduction 

For systems with long-range interaction the pair potential 
decays at large distances with a weak power law V(r) ~ 
r~ s , with s < d, the space dimension pQ. Examples are: 
gravity, whose statistical mechanics is made more com- 
plex by the unremovable singularity of the potential at 
the origin; Coulomb interactions, with the phenomenon 
of charge screening which allows the treatment of glob- 
ally neutral systems; dipolar media, that display the well 
known feature of shape dependence; vortices interacting 
with logarithmic potential in d = 2, for which Onsager 
first discussed microcanonical features, like the presence 
of negative temperatures, at the first Statphys meeting in 
Florence [2]. 

Long-range interactions can be made extensive, but 
are intrinsically non-additive. Let us consider the simplest 
model of ferromagnetic systems, the Curie- Weiss mean- 
field Hamiltonian 



H 



J 

2N 



(I) 



in which the spins Oi = ±1 are globally coupled with 
strength J > 0. The energy scales with system size, H ~ 
N, and an intensive energy density £ = Muin^^ H/N can 
be defined in the thermodynamic limit, but, due to the 



presence of iV 2 links, the sum of the energies of two sub- 
systems I and II is never giving the total energy Ej+u ^ 
Ei + En. Besides that, for mean-field models of this kind 
the set of accessible macrostates in the space of inten- 
sive parameters E and M — 2^£>i, can De non convex 
(see Fig.QJ. These two mathematical properties, that can 
be present only at finite iV for short-range interactions 
(s > d), have important physical consequences. The spe- 
cific heat [3] , and other quantities related to the curvature 
of the entropy, like susceptibility :4| , can become negative 
in certain energy ranges; jumps in temperature can be re- 
alized at first-order microcanonical transitions [5] ; broken 
ergodicity can be present both for finite N and in the ther- 
modynamic limit [5] (see Sect. [5]). All this is inscribed in- 
side the general framework of ensemble inequivalence [7|8j . 
We will briefly illustrate these features for a mean-field XY 
model with two-spin and four-spin interactions in Sect. [2j 

All the above is for equilibrium properties (i.e. max- 
imal entropy states), but systems with long-range inter- 
actions also show a very slow approach to equilibrium. 
For systems with short-range interactions it has been def- 
initely assessed that, appropriately selecting the initial 
state, the subsequent relaxation to equilibrium takes place 
on a finite time in the thermodynamic limit [9|. On the 
contrary, since the seminal paper of Lynden-Bell [3J, it 
has been proposed that systems with long-range interac- 
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tions can display a two-step relaxation. In a first stage, the 
system relaxes rapidly ( "violently" , according to Lynden- 
Bcll) to a quasi- stationary state whose lifetime increases 
with system size: two types of dependencies have been pro- 
posed, either power-law N s [TO] or logarithmic In N [B] , de- 
pending on some detailed property of the initial state. In a 
second stage, the system begins a slow approach to equilib- 
rium that may be either direct or proceed through succes- 
sive relaxations to different quasi-stationary states. The 
fact that the lifetime of quasi-stationary states diverges 
with system size allows to obtain a separation of the two 
time scales. This scenario, originally proposed for gravity, 
has been extended to the two-dimensional Euler equation, 
the original Onsager's system, by Chavanis A theo- 
retical proposal has been advanced by Lyndcn-Bcll in or- 
der to interpret the initial relaxation to a quasi-stationary 
state: that the system tries to maximize an entropy which 
takes into account additional "constraints" appearing in 
the thermodynamic limit, the most relevant of those be- 
ing the normalization of the one-body distribution func- 
tion [3] . Although the initial numerical tests of this theory 
gave some hope of success [T2], its predictions were sub- 
sequently disproved for several systems, even for the sim- 
ple one-dimensional self-gravitating system [13] and the 
theory fell into some discredit. We will show in Sect. [3] 
results for the so-called Hamiltonian Mean Field (HMF) 
model [14] (a mean-field XY model) for which Lynden- 
Bcll theory has a straightforward application and leads 
to predictions that are in reasonable agreement with nu- 
merical simulations [TB]. A recent domain of application 
of Lynden-Bell's ideas is to the free electron laser [15], 
where this approach allows to predict the features of the 
intensity saturation of the laser, simply knowing the ini- 
tial conditions, without explicitly solving the equations of 
motion, as was tipycally done before: in a sense, this con- 
stitutes the statistical mechanics of the free electron laser. 
This application is likely to produce in the future experi- 
mental tests of the features of quasi-stationary states, as 
we discuss in Sect. [6] 



Furthermore, the macrostates that maximize Lynden- 
Bell's entropy for a given energy depend on the initial 
magnetization of the HMF model. The system undergoes a 
nonequilibrium phase transition that can be of the second 
or the first order, and hence a nonequilibrium tricritical 
point is present (see Sect. 2]). It is interesting that the 
concepts of equilibrium phase transition can be extended 
to states that are not in equilibrium. It means that one 
can assume mixing limited to the phase-space visited by 
the trajectories on a time scale that is small with respect 
to the relaxation time to equilibrium. 



Sect. [7] is devoted to some final remarks and perspec- 
tives of applications to wave-particle systems and to dipo- 
lar media. 




M 



Fig. 1. The set of accessible macrostates in the (E,M) space 
can have a non-convex shape for systems with long-range in- 
teractions, such that if (Ei,M\) and (i?2,M2) can be realized 
macroscopically, this is not necessarily true for all the states 
joining these two along the straight dashed line. 



2 Phase diagram of a mean-field XY model 

Let us consider the following Hamiltonian 



Hxy — 



N n 

^ 2 



J 



N 
i=l 



K 



iV 



(E s *) 2 



(2) 



which can be thought as representing a system of N spins 
Si = (cos Oi, sin 6*j ) with all-to-all two-spin, J, and four- 
spin, K, interactions. A kinetic energy term is added, 
considering pi as canonically conjugate to Because of 
this addition, Hamiltonian ([2]) can also represent a sys- 
tem of unit mass particles moving on a circle without 
collisions, interacting only through a mean-field type po- 
tential. This model can be solved in both the canoni- 
cal and the microcanonical ensemble [T7]. The resulting 
phase diagram is shown in Fig. [2] For both ensembles a 
tricritical point is present, but its location is different in 
the two ensembles. The behavior of the order parame- 
ter m = limjv— >oo | ^iSil/N in the two ensembles is also 
shown, in order to highlight the striking difference in the 
predictions. The so-called caloric curve (kinetic tempera- 
ture vs. energy density £) is reported in Fig. [3] The micro- 
canonical ensemble (full line) predicts a region of negative 
specific heat, where kinetic temperature T decreases as 
the energy density is increased. Moreover, a temperature 
jump is present at the transition energy. 



3 Quasi-stationary states 

Hamiltonian ^ reduces, for K = and adding a constant 
to shift the energy of the ground state to zero, to that of 
the HMF model QJ] 



H 



HMF 



N r? 1 
9 9 N 



N 

TO- 



cos 



,)), (3) 



whose equilibrium properties are standard: the system un- 
dergoes a mean-field second order phase transition in both 
the microcanonical and the canonical ensemble at the en- 
ergy £ = 3/4, corresponding to the temperature T =1/2 
(see [IB] for a comprehensive recent review of the model) . 
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Fig. 2. Phase diagram of Hamiltonian @. The canonical 
second order transition line (full horizontal line starting at 
T/J — 1/2) becomes first order (dotted line) at the canonical 
tricritical point. The microcanonical second order transition 
line coincides with the canonical one up to K/J = 1/2 but 
it extends further towards the microcanonical tricritical point, 
located at K/J = 5/2. At this latter point, the transition line 
bifurcates in two first order microcanonical lines (dashed), cor- 
responding to a temperature jump. There are no microcanoni- 
cal macrostates for parameter values within the shaded region. 




Fig. 3. Caloric curve for K/J — 10. The full line is the theo- 
retical prediction in the microcanonical ensemble. The dashed 
line represents the first order phase transition in the canonical 
ensemble. The points are obtained from a molecular dynamics 
simulation of Hamiltonian ([2]) with TV = 100. 



In order to study nonequilibrium properties, it has been 
customary to prepare an initial state of the "water-bag" 
type where all the particles are uniformly distributed in 
a rectangular domain of width 2A9 and height 2Ap cen- 
tered around the origin in the single-particle phase space 
(9,p). Once the size of the domain is given, energy and 
magnetization are uniquely determined: mo = sin A9/A9, 
£ = (Ap) 2 /Q + (1 — m 2 ) )/2. During microcanonical time 
evolution, energy and total momentum J^. pi remain con- 
stant, but magnetization varies, and one expects that it 
reaches the equilibrium value compatible with the given 
energy. This is not what happens, as shown in Fig.[4j One 
observes an initial "violent" relaxation to a plateau value 
(closer and closer to m = as the number of particles 
is increased) , corresponding to the quasi-stationary state, 
followed by a second relaxation to equilibrium, which takes 
place on longer and longer times as TV increases. The life- 
time of the quasi-stationary state has been fitted with a 
power-law TV 1 - 7 in Ref. [TU] , 

How does one explain such a behaviour? It is crucial to 
realize that the mean-field dynamics ([3|) is well represented 
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Fig. 4. Magnetization m vs. time t for the HMF model with en- 
ergy £ = 0.69, vanishing total momentum and "water-bag" ini- 
tial condition with mo = (AO = n). The value of TV increases 
from left to right: TV = 10 2 ,10 3 ,2 x 10 3 ,5 x 10 3 , 10 4 ,2 x 10 4 . 



in the large TV limit by the following Vlasov equation 



Bt P d9 d9 dp 



(4) 



where f(9 1 p,t) is the single-particle distribution function 
and the potential V(9) is given by 

V(6) [f] = 1 - m x [f] cos(9) - m y [f] sm(9), 

m x [f] = f f(0,p,t) cos9d9dp, 



[/] 



f(9,p,t) sm9d9dp. 



Besides energy, Vlasov equation also conserves the norm 
of the distribution function J f(9,p,t)d9dp = 1, which 
for the "water-bag" initial condition we have considered 
implies that the distribution function remains two-level 
(0,/ = l/(AA9Ap)) at all times. To take this into ac- 
count, Lyndcn-Bcll has proposed [3] that the system adapts 
itself on macrostates that maximize the following "fermionic" 
entropy 



slb{J) 



dpd9 



fo /o 



1-^ 111) 

Jo 



(5) 

where / is the coarse-grained distribution function. The 
Pauli principle that is implicit in this approach refers to 
the fact that a "fluid element" in the (9,p) /Lt-space can- 
not occupy a cell which is already occupied by another 
fluid element, just because the distribution must remain 
two-level in the course of time. It turns out that the max- 
imization of SLs(f) at fixed energy, momentum and norm 
can be explicitly performed for the HMF model, giving 
the following solution 



fQSs{0,p) 



.fo 



e /3(p 2 /2-m[/ QSS ]cos0)+Ap+a _|_ y 



(G) 



where (3, A and a are Lagrange multipliers correspond- 
ing to the conservation of energy, momentum and norm. 
As it is clear from equation ©, the quasi-stationary dis- 
tribution has to be determined self-consistently, since it 
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Fig. 5. (color online) Comparison of the theoretical predic- 
tions for the quasi-stationary momentum distribution (dashed 
line) with the numerical simulations (points) performed using 
Hamiltonian © for £ = 0.69 and mo = 0.3 (a), mo = 0.5 (b), 
mo = 0.7 (c) . The case mo = 0.3 is also shown in lin-lin scale 
in panel (d). For all these cases m[fQss] = (homogeneous 
quasi-stationary states) . 



depends on the distribution itself through the magneti- 
zation to[/qss]. We have chosen the subscript QSS to 
mean quasi-stationary-state, because we propose [161 that 
the states that maximize Lynden-Bell entropy are indeed 
the quasi-stationary states observed, e.g., in the numeri- 
cal simulation reported in Fig. [4] The magnetization in the 
QSS, m[fQss] — niQss, and the values of the Lagrange 
multipliers are obtained by solving numerically a set of 
implicit equations [16j . once the energy £ and the initial 
magnetization mo are given (or alternatively (A9,Ap)). 
Since we look for solutions whose total momentum van- 
ishes, A = 0. In Fig. [5] we show the comparison of the 
predictions of the theory concerning momentum distribu- 
tions with the numerical simulations performed integrat- 
ing the equations of motion given by Hamiltonian (|3|) with 
iV = 10000. The agreement is quite good, if one takes into 
account that the predictions have no adjustable parame- 
ter and are strictly determined by the choice of the ini- 
tial condition. However, expecially the plot in lin-lin scale 
(panel (d) in Fig. [5} reveals that our theory is unable to 
reproduce the double bump obtained in numerical simu- 
lations. The presence of this bumpy feature is confirmed 
by simulating the Vlasov equation |4]) (see Fig. [6]) [20] . 
The double bump in momentum distribution is the result 
of the presence of two sliding resonances (the two vortices 
in the lower right panel), for which an explanation that 
takes into account the specific dynamical properties of the 
model is necessary |16l20j . 

Quasi-stationary states have been demonstrated to be 
robust to both the application of and external bath [21] 
and to the addition of a small nearest neighbor interac- 
tion [22 . The Vlasov equation is exact in the N — > oo 
limit, and people have tried to address the question of fi- 
nite N corrections. In this respects the most interesting 
progress is registered in papers by Bouchet and Dauxois 
[23] and Chavanis [24]. For the sake of commpleteness, 
it should be mentioned that a different approach, based 
on Tsallis statistics, has been proposed to describe quasi- 



Fig. 6. (color online) Time evolution of the single-particle dis- 
tribution function according to the Vlasov equation ([4]) for 
a rectangular "water-bag" initial state with £ — 0.69 and 
m = 0.5. 
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Fig. 7. (color online) Theoretical phase diagram on the con- 
trol parameter plane (mo,£): second order phase transition 
line (dashed); first order phase transition line (full); tricritical 
point (full dot). Inset: magnification of the first order phase 
transition region and limits of the metastability region (dash- 
dotted). 



stationary states [25] ■ Recently, these authors have con- 
centrated their attention on the behaviour of correlations 
in the single particle dynamics, that can give rise to non 
Gaussian distributions for sums of variables (generalized 
central limit theorems). 



4 Nonequilibrium tricritical point 

The maximization of Lynden-Bell entropy reserves an- 
other surprise. Since this variational problem introduces 
another control parameter besides energy, the magnetiza- 
tion mo of the initial state, one obtains a phase transition 
at an energy that depends on mo. The transition energy 
coincides with that given by Boltzmann entropy (to which 
Lynden-Bell's entropy reduces in the diluted limit / <C /o) 
only for mo = 1. Not only, below mo = 0.61 the transition 
becomes first order. The full phase diagram is plotted in 
Fig. [7] A nonequilibrium tricritical point is present in the 
phase diagram. The specific heat is negative along all the 
transition line and a temperature jump appears when the 
transition change to first order. 

Numerical simulations confirm this extremely inter- 
esting theoretical finding [26] . For instance, we plot in 
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Fig. 8. (color online) rriQss as a function of £ for mo = 0.30 
where the phase transition is second order. The full line is the 
theoretical prediction and the points are simulations of the 
HMF model for which the number of initial realizations is 10 
(N = 10 3 ), 10 4 (N = 10 4 ), 10 2 (N = 10 6 ) and the averaging 
time 20 < t < 100. 
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Fig. 9. (color online) rriQss as a function of £ for mo = 0.05 
where the phase transition is first order. The results are plotted 
as in Fig. [8] 

Figs. I8|9I the order parameter itiqss as a function of en- 
ergy density at fixed mo, finding signature of a second 
order (Fig. [S} or first order (Fig. [9|) phase transition. As 
previously explained, these are "short-time" phase transi- 
tions, at variance with the "long-time" phase transitions 
of equilibrium statistical mechanics. Hence, in order to re- 
veal them numerically, one has to average over a short 
initial time and take many initial instances. 



5 Broken ergodicity 

After this excursion in nonequilibrium, let us come back to 
equilibrium properties. The XY model ([2]), for some spe- 
cific choice of antiferromagnetic J < and ferromagnetic 
K > couplings has energy density values where regions 
with different magnetization m are disconnected if energy 
is fixed (as shown in Fig. [T]) [57]. An example is given in 
Fig. [lOl where we show two different situations. In panel 
a) the "standard" ergodicity breaking associated with a 
phase transition is revealead by successive magnetization 
switches between local entropy maxima: the number of 
particles is small (N = 20) such that the entropy barrier 
(see inset) is not insurmontable. In panel b) we display 
the new type of microcanonical broken ergodicity first dis- 
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t 



Fig. 10. Time evolution of the magnetization m of model ([2]) 
for: a) £ = 0.1, J = -1 and K = 8; b) £ = 0.1077, J = -1 and 
K = 3, obtained integrating the equations of motion derived 
from Hamiltonian (J2j) with TV = 20. In panel a) magnetization 
switches between the two most probable values (see the inset 
for the dependence of entropy on m) . In panel b) , two different 
initial conditions are plotted simultaneously, corresponding to 
two different values of the initial magnetization, mo = 0.1 and 
mo = 0.98. The inset shows the entropy, which now vanishes 
in the interval [m_,m+]. Here m_ ~ 0.192 and m+ ~ 0.794, 
which are shown to bound from above (dotted line) and from 
below (dashed line) the magnetization of the two initial condi- 
tions. 

cussed in Refs. [6,28]: now magnetization cannot switch 
because there is a gap [m_,m. + ] with no macrostates. 



6 Application to the free electron laser 

Free-Electron Lasers (FELs) are coherent and tunable ra- 
diation sources, which differ from conventional lasers be- 
cause they use a relativistic electron beam as their las- 
ing medium, hence the term "free-electron" . The physi- 
cal mechanism responsible for light emission and ampli- 
fication is the interaction between the relativistic elec- 
tron beam, a magnetostatic periodic field generated in 
the undulator and an optical wave co-propagating with 
the electrons. Due to the presence of the magnetic field, 
electrons are forced to follow oscillating (tipycally sinu- 
soidal) trajectories and emit synchrotron radiation. This 
spontaneous emission is then amplified along the undula- 
tor until the laser effect is reached. Among the different 
schemes, single-pass high-gain FELs are currently attract- 
ing a growing interest. Basic features of the system dy- 
namics are successfully captured by a simple Hamiltonian 
model introduced by Colson and Bonifacio [2"5] 
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dA 1 



(7) 



where z = 2k w pzj 2 / (7)0 is the rescaled longitudinal coor- 
dinate, which plays the role of time. Here, p = 7,~ 1 (xrr- £ ) 2 ' 3 



similarities in the behaviours of the two models. Indeed, 
similarly to the HMF model, Hamiltonian presents a 
standard second order phase transition of the mean field 
type. As for the HMF, FEL dynamics is well represented 
in the N — > 00 limit by the following Vlasov equation 



is the so-called Pierce parameter, (7)0 the mean energy of 
the electrons at the undulator's entrance, k w = 2ir/\ w the 



wave number of the undulator, uj v 



(4ire 2 n/m) 1 / 2 the 
plasma frequency, c the speed of light, n the total elec- 
tron number density, e and m respectively the charge and 
mass of one electron. Furthermore, a w = eB w /(k w mc : ), 

where B w is the rms peak undulator field. Here "f r = 
1/2 

(\ w (l + a 2 u )/2X) is the resonant energy, X w and A be- 
ing respectively the period of the undulator and the wave- 
lenght of the radiation field. Introducing the wavenum- 
ber k of the FEL radiation, the two canonically conju- 
gated variables are (9,p), denned as 9 = (k + k w )z — 
2Spk w z-f 2 /('j)l andp = (7- (7>o)/(p(7>o)- corresponds 
to the phase of the electrons with respect to the pon- 
deromotive wave. The complex amplitude A = A x + iA y 
represents the scaled field, transversal to z. Finally, the 
detuning parameter is given by S = ((7)0 — 7r)/(2p7r)> 
and measures the average relative deviation from the res- 
onance condition. This succinct, but detailed, description 
of the model should transmit the idea that, by adjusting 
model parameters, one can really simulate realistic exper- 
imental situations. 

It can be easily checked that model © can be derived 
from the Hamiltonian 



H 



FEL 



N 2 

Y*L-N5I- 



A' 



V), (8) 



where N is the number of electrons and the intensity / and 
the phase cp of the wave are related to A = A. x + iA y 
\ I e~ tip . In addition to the energy, the total momentum 
P = Y^j Pj + AAA* is also a conserved quantity. 

It should also be mentioned that, remarkably, this sim- 
plified formulation applies to other physical systems, pro- 
vided an identification of the variables involved is per- 
formed. As an example, consider the electron beam-plasma 
instability. When a weak electron beam is injected into a 
thermal plasma, electrostatic modes at the plasma fre- 
quency (Langmuir modes) are destabilized. The interac- 
tion of the Langmuir waves and the electrons constitut- 
ing the beam can be studied in the framework of a self- 
consistent Hamiltonian picture [31] , formally equivalent to 
the one in [29] . In a recent paper [32] we have established 
a bridge between these two areas of investigation and ex- 
ploited the connection to derive a reduced Hamiltonian 
model to characterize the saturated dynamics of the laser. 

There are many similarities between model (JSJ and the 
HMF model ([3]), and indeed the canonical free energy of 
one model can be exactly mapped onto the other [5U] , 
While in the HMF model particles interact directly, in the 
FEL dynamics electrons interact only through the field 
A, whose dynamics depend in turn on those of the elec- 
trons. Therefore, already on this basis one could expect 



df 

dz 
dA x 
dz 

dAy 

dz 



df 



= - 5A v + 



2{A X cos 6 — A y sin( 



1 



) — 

dp' 



SA, 



1 

2^ 



/ cos 9 d9dp, 
f sin 9 d9dp . 



(9) 



Hence, one can similarly expect the existence of quasi- 
stationary states 15J. A direct simulation of Eqs. is 
shown in Fig. [TT] for N — 10 4 electrons. Initially, the elec- 
trons are uniformly distributed in the 9 interval [— 7r,7r] 
and their momentum is also uniformly spread in [~Ap, Ap] 
in a "water-bag" distribution: this initial state has a phys- 
ical meaning, because it is very close to the state in which, 
experimentally, the electrons are injected into the undu- 
lator. The laser intensity / is initially set to a small value, 
in order to switch on the instability, whose initial growth 
is exponential. On a short time, the intensity reaches a 
saturation level and performs wide oscillations around it, 
but these oscillations dump with time, and a well de- 
fined asymptotic intensity is reached. This process is the 
equivalent of Lynden-Bell's "violent relaxation" . The N- 
dependence of the saturation level is shown in the inset, 
where three different values of TV are considered. For the 
smaller value, N = 100 (curve 3), the saturated inten- 
sity is larger than for higher N values. When N increases 
the system remains trapped in the quasi-stationary state 
for a longer and longer time. The second relaxation to 
the larger intensity value corresponds to the relaxation to 
Boltzmann-Gibbs equilibrium. This latter relaxation has 
no experimental relevance for free electron lasers, because 
it would take place for enormous undulator lengths. A 
similar behavior would be observed for the "bunching pa- 
rameter" b — I J2 n ex P(i@n)\/N, which is the equivalent of 
the magnetization m of the HMF model. 

The distribution fQss{6, p) in the quasi-stationary state 
can be obtained by maximizing Lyndcn-Bcll entropy (JSJ), 
keeping energy and momentum P fixed. The result of this 
maximization procedure is very similar to that of the HMF 



Jqss(9,p) 



fo 



s (3(p 2 /2+2A Q ss sin0)+\p+a _|_ ^ ' 



(10) 



with 



Aqss = \J Iqss — 



/36 -A 



dpd9 sm9f QSS (9,p), (11) 



and /3, A and a are the usual Lagrange multipliers, whose 
value is determined by the initial condition. 

The dependence of both the intensity and the bunch- 
ing parameter in the quasi-stationary state on the de- 
tuning parameter S for an initially homogeneous state 
(A9 = 7r) with zero momentum dispersion Ap = is 
shown in Fig. (|12[). 
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20 40 60 80 100 120 140 160 180 200 
2 

Fig. 11. Typical evolution of the radiation intensity of a free 
electron laser using Eqs. J7J; the detuning 5 is set to 0, the 
energy per electron H/N = 0.2 and TV = 10 4 electrons are 
simulated. The inset presents averaged simulations on longer 
times for different values of TV: 5 ■ 10 3 (curve 1), 400 (curve 2) 
and 100 (curve 3). 
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Fig. 12. Laser intensity Iqss and bunching parameter bQss 
in the quasi-stationary state as a function of detuning 5: the- 
oretical prediction (solid and dashed lines, respectively) and 
simulations (symbols). The dotted vertical line, 8 — S c ~ 1.9, 
represents the transition from the low to the high-gain regime. 



If FEL's Hamiltonian is equivalent to HMF, one might 
wonder if situations exist which produce more complex 
phase diagrams and dynamical evolutions. Indeed, in the 
regime where the beam current and the emittance is small, 
FEL's dynamics is better described by a model where 
many planar waves interact with the beam 

604 



ihO. 



A* h e- ih9 ') 



— ihOj 



dz 



(12) 



where the complex amplitudes Ah = + iA v h represent 
the scaled field, transversal to z, and the coupling pa- 




Fig. 13. (color online) Phase diagram in the (F 1 ,F^) parame- 
ter space for the two- wave model (|12p . corresponding to a FEL 
with an initial cold homogeneous beam. Lynden-Bell's theory 
predicts a transition from a region Z 1 (shaded) dominated by 
the first wave to one, Z3, dominated by its odd harmonic. 
Above the dashed line the short time growth rate of | A3 1 is 
predicted to be larger than that of |Ai| by a linear theory 



rameters Fh depend on the experimental setup. We have 
analysed in some detail [33j only the case with two odd 
harmonics h = 1,3, performing the same analysis of the 
quasi-statonary states as for Eqs. 0. The initial condition 
is again a homogeneous almost monocromatic beam, but 
now there are two "order parameters" , |Ai| and | A3 1 , with 
the corresponding bunchings b\ = | J^ . exp(— i6j)/N and 
b 3 = I £\ eKp(-i30j)/N. In region Z x (shaded) of Fig. QJ 
the quasi-stationary state is dominated by the first wave: 
hence |Ai| > and | A3 1 = 0. On the contrary, in region 
Z3, I A3 1 > and |Ai| = 0. The full line corresponds to a 
nonequilibrium phase transition of the first order. 

However, a linear analysis show that in the region be- 
tween the dashed line and the transition line to Z 3 (inside 
the Zx region) the growth rate of |A 3 | is larger than that 
of I Ai I = 0. Therefore, at short time the growth of | A3 1 is 
expected to be faster and the initial evolution of the FEL 
is dominated by the first odd harmonic. However, this is 
not the maximal Lynden-Bell's entropy state, and on a 
later time, the system is finally going to be dominated by 
the first wave. This is indeed what happens, as shown in 
Fig. M 



7 Conclusions and perspectives 

We have discussed both equilibrium and nonequilibrium 
properties of systems with long range interactions with ref- 
erence to mean-field models. Most of these features should 
extend to cases where forces weakly decay |H30j . but this 
problem remains to be seriously investigated. In this per- 
spective, the general classification of phase transitions ob- 
tained in Ref. [34] could be extremely useful. 

The domain of application of these ideas is vast and the 
perspective to perform key experiments in the near future 
is realistic. We have discussed in some detail in Sect. [6] a 
realistic model of the free electron laser. However, we can 
foresee applications to all systems where particles interact 
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Fig. 14. (color online) Upper panel: Evolution of the inten- 
sities | Ai| 2 and | A3 1 2 showing an initial growth of the first 
odd harmonic followed by the final relaxation to the maximal 
Lynden-Bell's entropy state dominated by the first wave. Lower 
panel: time evolution of Lynden-Bell's entropy 

with waves in a self-consistent way. As an example, let 
us mention collective atomic recoil lasers, where ultracold 
atoms interact with the optical modes of the cavity, deter- 
mining collective effects and first and second order phase 
transitions |35j . 

Another domain where experimental applications are 
envisageable is that of layered spin structures, where dipo- 
lar interactions dominate over Heisenberg exchange |36j . 
In this case, a mean-field term of the Hamiltonian is shown 
to depend on the shape of the sample, determining the 
presence or absence of phase transitions [37] . 

Most of my knowledge of the physics of systems with long- 
range interactions has been shaped by J. Barre, F. Bouchet, 
P.H. Chavanis, T. Dauxois, D. Fanelli D.H.E. Gross, D. Mukamel, 
Y.Y. Yamaguchi. I also thank all the coauthors of my papers 
on this subject, with whom I have entertained endless discus- 
sions. This work is funded by the PRIN05 grant Dynamics and 
thermodynamics of systems with long-range interactions. 
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